REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this 
collection  of  information,  including  suggestions  for  reducing  this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and 
Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person 
shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR 
FORM  TO  THE  ABOVE  ADDRESS. 


2.  REPORT  TYPE 

Journal  Article 


1.  REPORT  DATE  (DD-MM-YYYY) 
September  2013 


4.  TITLE  AND  SUBTITLE 

Helium  Nanodroplet  Isolation  and  Infrared  Spectroscopy  of  the  Isolated  Ion- 
Pair  1  -ethyl-3  -methy  limidazolium  bis(trifluoromethylsulfonyl)imide 


6.  AUTHOR(S) 

E.I.  Obi,  C.M.  Leavitt,  P.L.Raston,  C.P.  Moradi,  S.D  Flynn,  J.A.  Boatz,  S.D. 
Chambreau,  G.E.  Douberly 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Research  Laboratory  (AFMC) 

AFRL/RQRP 
10  E.  Saturn  Blvd. 

Edwards  AFB  CA  93524-7680 


9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Research  Laboratory  (AFMC) 

AFRL/RQR 
5  Pollux  Drive 

Edwards  AFB  CA  93524-7048 


12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Distribution  A:  Approved  for  Public  Release;  Distribution  Unlimited.  PA#  13471 


3.  DATES  COVERED  (From  -  To) 
September  2013-  October  2013 


5a.  CONTRACT  NUMBER 

In-House 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 

Q0RA 


8.  PERFORMING  ORGANIZATION 
REPORT  NO. 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 


11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

AFRL-RQ-ED- J  A-2 013-208 


13.  SUPPLEMENTARY  NOTES 

Journal  Article  submitted  to  Journal  of  Physical  Chemistry. 


14.  ABSTRACT 

The  ionic  liquid  l-ethyl-3-methylimidazolium  bis(trifluoromethylsulfonyl)imide  is  vaporized  at  420  K,  and  the  ion-pair  constituents  are 
entrained  in  a  beam  of  liquid  He  nanodroplets  and  cooled  to  0.4  K.  The  vapor  pressure  is  optimized  such  that  each  He  droplet  picks  up  a 
single  ion-pair  from  the  gas  phase.  Infrared  spectroscopy  in  the  CH  stretch  region  reveals  bands  that  are  assigned  to  intact  ion-pairs  on  the 
basis  of  comparisons  to  ab  initio  harmonic  frequency  computations  of  23  low  energy  isomers.  The  He  droplet  spectrum  is  consistent  with  a 
weighted  sum  of  the  computed  harmonic  spectra,  in  which  the  weights  are  detennined  from  ab  initio  computations  of  the  relative  free 
energies  at  420  K.  Anharmonic  resonance  polyads  in  the  CH  stretch  region  are  treated  explicitly,  which  improves  the  agreement  between  the 
experiment  and  computed  spectra  for  ion-pairs.  For  isomers  having  a  strong  cation... anion  hydrogen  bonding  interaction,  the  imidazolium 
C(2)-H  stretch  fundamental  is  shifted  to  lower  energy  and  into  resonance  with  the  overtones  and  combination  of  the  imidazolium  ring 
stretching  modes,  resulting  in  a  spectral  complexity  in  the  CH  stretch  region  that  is  fully  resolved  in  the  He  droplet  spectrum.  The  assignment 
of  the  infrared  spectrum  to  ion-pairs  is  confirmed  through  polarization  spectroscopy  measurements  that  reveal  the  permanent  electric  dipole 
moment  of  the  He-solvated  species  to  be  1 1  2D.  The  computed  permanent  electric  dipole  moments  for  the  low  energy  isomers  of  the 
[emim+][Tf2N-]  ion  pairs  fall  in  the  range  9-13  D,  whereas  the  computed  dipole  moments  of  decomposition  products  of  the  ionic  liquid  are 
less  than  4.3  D. 


15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 


a.  REPORT 


Unclassified 


b.  ABSTRACT  c.  THIS  PAGE 


Unclassified  Unclassified 


17.  LIMITATION 

OF  ABSTRACT 

18.  NUMBER 

OF  PAGES 

SAR 

31 

19a.  NAME  OF 
RESPONSIBLE  PERSON 

Tom  Hawkins 


19b.  TELEPHONE  NO 

(include  area  code) 

661-525-5449 

Standard  Form 
298  (Rev.  8-98) 

Prescribed  by  ANSI 


Helium  Nanodroplet  Isolation  and  Infrared  Spectroscopy  of  the  Isolated  Ion- 
Pair  l-ethyl-3-methylimidazolium  bis(trifluoromethylsulfonyl)imide 


Emmanuel  I.  Obi',  Christopher  M.  Leavitt^,  Paul  L.  Raston',  Christopher  P.  Moradi',  Steven  D. 
Flynn’,  Jerry  A.  Boatz^’*,  Steven  D.  ChambrcarL’*  and  Gary  E.  Doubcrly1* 


Department  of  Chemistry,  University  of  Georgia,  Athens,  GA  30602 
§Air  Force  Research  Laboratory,  Edwards  Air  Force  Base,  CA  93524 
ARC  Incorporated,  Edwards  Air  Force  Base,  CA  93524 


Corresponding  authors:  jerry.boatz@us.af.mil  steven, chambreau. ctr@us.af.mil +'* 

douberly@iiga.edu 


Abstract 

The  ionic  liquid  1 -ethyl-3 -methylimidazolium  bis(trifluoromethylsulfonyl)imide  is  vaporized  at  420  K, 
and  the  ion-pair  constituents  are  entrained  in  a  beam  of  liquid  He  nanodroplets  and  cooled  to  0.4  K.  The 
vapor  pressure  is  optimized  such  that  each  He  droplet  picks  up  a  single  ion-pair  from  the  gas  phase. 
Infrared  spectroscopy  in  the  CH  stretch  region  reveals  bands  that  are  assigned  to  intact  ion-pairs  on  the 
basis  of  comparisons  to  ab  initio  harmonic  frequency  computations  of  23  low  energy  isomers.  The  He 
droplet  spectrum  is  consistent  with  a  weighted  sum  of  the  computed  harmonic  spectra,  in  which  the 
weights  are  determined  from  ab  initio  computations  of  the  relative  free  energies  at  420  K.  Anharmonic 
resonance  polyads  in  the  CH  stretch  region  are  treated  explicitly,  which  improves  the  agreement  between 
the  experiment  and  computed  spectra  for  ion-pairs.  For  isomers  having  a  strong  cation  anion  hydrogen 
bonding  interaction,  the  imidazolium  C(2)-H  stretch  fundamental  is  shifted  to  lower  energy  and  into 
resonance  with  the  overtones  and  combination  of  the  imidazolium  ring  stretching  modes,  resulting  in  a 
spectral  complexity  in  the  CH  stretch  region  that  is  fully  resolved  in  the  He  droplet  spectrum.  The 
assignment  of  the  infrared  spectrum  to  ion-pairs  is  confirmed  through  polarization  spectroscopy 
measurements  that  reveal  the  permanent  electric  dipole  moment  of  the  He-solvated  species  to  be  1 1±2  D. 
The  computed  permanent  electric  dipole  moments  for  the  low  energy  isomers  of  the  [emim+][Tf2N']  ion- 
pairs  fall  in  the  range  9-13  D,  whereas  the  computed  dipole  moments  of  decomposition  products  of  the 
ionic  liquid  are  less  than  4.3  D. 
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Introduction 


Ionic  liquids  are  low  melting  point  salts  (<  373  K)  whose  bulky  constituents  lead  to 
characteristically  long  cation-anion  distances  and  low  lattice  energies,  in  comparison  to  common 
ionic  salts,  such  as  NaCl.  ’  For  many  cation-anion  combinations,  ionic  liquids  possess 
properties  such  as  low  volatility,  thermal  stability,  and  a  marked  stability  towards  air  and  water. 
In  principle,  these  physical  and  chemical  properties  can  be  tuned  by  the  choice  of  cation-anion 
combination,  making  ionic  liquids  potentially  useful  for  batteries,  '  solar  cells,  ’  fuel  cells,  ’ 
and  as  “green”  solvents10  for  organic  synthesis  and  chromatography.11  Ionic  liquids  have  also 
been  shown  to  exhibit  considerable  promise  as  fuels  for  electrospray  ion  thrusters  in  space 
propulsion  applications.12'14  Chambreau  and  co-workers  reported  a  class  of  room  temperature 
ionic  liquids  (RTILs)  that  spontaneously  ignite  upon  the  introduction  of  a  strong  oxidizer,  such 
as  fuming  nitric  acid,15, 16  and  other  groups  have  more  recently  reported  the  ignition  of  RTILs 
with  nitric  acid,  '  nitrogen  tetroxide,"  and  hydrogen  peroxide.'  This  observed  hypergolic 
activity  suggests  that  a  properly  chosen  RTIL  /  oxidizer  combination  could  replace  the  highly 
volatile  and  environmentally  toxic  hydrazine-based  hypergolic  fuels  currently  used  in 
bipropellant  rocket  engines.16  Given  the  vast  number  of  cation-anion  combinations,  it  is 
impractical  to  experimentally  characterize  even  a  small  fraction  of  all  possible  RTILs.  The 
optimization  of  these  systems  for  use  in  the  above  applications  will  therefore  require  predictive 
models  of  RTIL  properties.  The  molecular  level  understanding  on  which  these  models  must  be 
based  will  be  partially  derived  from  experimental  methods  capable  of  probing  the  fundamental 
aspects  of  RTILs  isolated  in  vacuo,  such  as  high  resolution  spectroscopy.  We  report  here  the 
vibrational  spectrum  and  dipole  moment  of  a  prototype  RTIL  that  is  vaporized  and  entrained  in  a 
beam  of  liquid  He  nanodroplets. 

Until  recently,  it  was  presumed  that  RTILs  possess  essentially  no  vapor  pressure,  due  to 
the  ionic  nature  of  the  liquid.  Earle  and  co-workers  were  the  first  to  demonstrate,  however,  that 
the  vacuum  distillation  of  certain  families  of  aprotic  RTILs  could  be  achieved  below  the 
decomposition  temperature  of  the  material."  This  result  pointed  towards  the  possibility  of  gas- 
phase  experiments,  and  many  reports  have  since  been  presented  that  investigate  both  the  nature 
of  the  ionic  liquid  vapor  and  the  vaporization  mechanism.  For  many  of  these  gas  phase  studies, 
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the  ionic  liquid  [1  -ethyl-3 -me thyhmidazolium][bis(trifluoromethylsulfonyl)imide]  (abbreviated 
[emim+][Tf2N'];  see  Scheme  1  and  Figure  1)  was  used  due  to  its  low  viscosity  and  high  thennal 
stability  at  the  temperatures  required  to  produce  a  sufficient  vapor  pressure.  This  prototypical 

23 

RTIL  has  been  studied  in  vacuo  with  line  of  sight  electron  impact  mass  spectrometry  (MS),' 
VUV  photoionization  MS,24  25  FTICR-MS,26  selected-ion  flow  tube  MS.27  These  indirect  probes 
of  the  neutral  gas-phase  constituents  are  all  consistent  with  the  vaporization  of  individual  ion- 
pairs  with  little  or  no  decomposition.  Consistent  with  this,  enthalpies  of  vaporization  (AHvap) 
were  computed  to  be  lowest  for  ion-pairs  ,  and  computations  of  AHvap  '  for  [emim  ][Tf2N"] 

30 

agree  well  with  experiment.  Furthermore,  the  valence  binding  energy  spectrum  obtained  from 

9Q 

photoelectron  spectroscopy  and  the  photoionization  appearance  energy  measured  with  tunable 
VUV  synchrotron  radiation  are  in  agreement  with  computed  valence  binding  energies  and 
adiabatic  ionization  potentials  of  individual  ion-pairs.  Collectively,  these  studies  of  vaporized 
1  -alkyl-3 -methylimidizolium  based  ionic  liquids  provide  strong  evidence  for  intact  ion-pairs.32 
Two  studies  have  directly  detected  the  intact  ion-pairs  of  a  similar  dialkylimidazolium-based 
RTIL."’  Nevertheless,  there  has  been  a  lack  of  experimental  studies  capable  of  probing  the 
detailed  geometric  structure  of  these  gas-phase,  isolated  ion-pairs. 


SCHEME  1 
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There  is  a  growing  literature  of  ab  initio  and  DFT  computations  of  the  geometric 
structure  and  electronic  properties  of  gas-phase,  ionic  liquid  ion-pairs.35’43  Computations  predict 
several  conformational  isomers  associated  with  isolated  [emim+][Tf2N‘]  ion-pairs,42’ 44  and  even 
for  this  highly  studied  system,  there  still  remains  significant  debate  over  the  preferred  geometric 
structure.  There  is  generally  little  agreement  on  the  relative  role  of  charge-charge,  hydrogen 
bonding  and  pi-stacking  interactions  in  determining  the  structure  of  isolated  ion-pairs.45  Tsuzuki 
and  co-workers  have  downplayed  the  importance  of  the  orientation  of  the  [TffNT]  anion  relative 
to  the  acidic  imidazolium  C(2)H  bond  and  suggested  instead  that  the  dominant  charge-charge 
interactions  result  in  a  rather  flat  potential  surface  with  many  isomers  having  similar  binding 

energies.  On  the  other  hand,  hydrogen  bonded  structures  with  multiple  C-FT  X  hydrogen 
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bonds  (X  =  N,  F,  or  O)  are  often  invoked  in  order  to  assign  vibrational  bands  in  condensed  phase 
IR  spectroscopy  experiments  of  [emim+] [T fiN"] . 46-50  For  example,  in  cryogenic  neon  matrix 
isolation  experiments,44'51'52  it  was  concluded  that  [emim+][Tf2N‘]  structures  having  strong  in 
plane  C(2)H+'N"  hydrogen  bonds  were  present  in  large  abundance  on  the  basis  of  comparisons 
between  the  experimental  FTIR  spectra  in  the  52  to  1700  cm"1  region  and  ab  initio  harmonic 
frequency  computations.44 

The  analysis  of  many  of  these  condensed  phase  studies  (liquid  room  temperature  or 
cryogenic  matrix  isolation)  have  relied  only  on  comparisons  between  FTIR  spectra  and  DFT  or 
ab  initio  (MP2  level  of  theory)  harmonic  frequency  computations.  Interestingly,  a  recent 
computation  of  the  cubic  force  constants  for  the  bare  [emim4  ]  cation  suggests  the  presence  of 
strong  anharmonic  resonances  in  the  CH  stretch  region/  This  motivated  Lassegues  and  co¬ 
workers  to  reassign  the  CH  stretch  region  of  condensed  phase  FTIR  spectra  to  a  class  of 
structures  that  does  not  involve  strong,  directional  C(2)H  '  anion  hydrogen  bonds.53' 54  The 
conclusions  drawn  by  these  authors  from  these  anharmonic  frequency  computations  contradict 
those  of  several  previous  reports  that  relied  on  the  relative  intensities  of  bands  in  the  CH  stretch 
region  to  make  inferences  about  the  structural  properties  of  the  liquid.46  Given  this 
disagreement55' 56  that  persists  in  the  literature,  there  is  a  clear  need  for  high  resolution  IR 
spectroscopic  methods  that  probe  the  structure  of  isolated  gas  phase  ion-pairs.  Helium 
Nanodroplet  Isolation  is  a  well  suited  method  for  this  purpose.  "  The  large  capture  cross 
sections  associated  with  He  droplets  allows  for  the  pick-up,  trapping,  cooling,  and  spectroscopic 
interrogation  of  fragile  species  that  have  relatively  low  vapor  pressures.  Over  a  1  cm  path 
length,  only  «10' 1  molecules/cm3  (10"5  Torr)  are  required  to  dope  the  droplets  with  single  RTIL 
constituent  molecules.  The  weakly  perturbing  He  solvent  allows  for  highly  resolved  vibrational 
spectroscopy  studies  of  the  solvated  molecular  species,  which  can  be  compared  directly  to  the 
predictions  of  quantum  chemistry.  Indeed,  when  comparisons  are  available,  the  vibrational  band 
origins  of  He-solvated  molecules  and  molecular  complexes  differ  little  from  those  measured  in 
the  gas  phase  (~1  cm"1  or  less).57 

As  discussed  in  the  present  report,  the  prototypical  RTIL,  [emim+][Tf2N"],  was  vaporized 
at  410  K  and  entrained  in  liquid  He  droplets.  On  the  basis  of  previous  reports  employing  this 
methodology,63' 64  we  expect  the  gas  phase  thennal  distribution  of  [emim+][Tf2N"]  isomers  to  be 
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“frozen-out”  due  to  the  rapid  dissipation  of  the  molecular  internal  energy  upon  He  solvation. 

The  IR  spectrum  in  the  CH  stretch  region  is  consistent  with  a  420  K  thermal  distribution  of  intact 
ion-pairs,  having  both  hydrogen  bonded  and  stacked  structures.  The  high  resolution  associated 
with  this  measurement  allows  for  a  direct  comparison  of  the  observed  spectrum  to  ab  initio 
frequency  computations  of  23  [emim+][Tf2N‘]  isomers.  For  satisfactory  agreement  between 
experiment  and  theory,  the  anharmonic  coupling  between  ring  stretching  overtones  and  CH 
stretching  fundamentals  must  be  explicitly  considered,  indicating  that  strong  anharmonic 

53 

resonances  persist  in  this  spectral  region,  as  originally  suggested  by  Lassegues  and  co-workers. 
Furthennore,  the  intensity  of  a  well-resolved  IR  band  was  measured  as  the  He-solvated  molecule 
was  oriented  in  a  lab-frame  electric  field.  This  polarization  spectroscopy  measurement  is 
consistent  with  an  [emim+][Tf2N‘]  pennanent  electric  dipole  moment  of  1 1±2  D,  which 
corroborates  the  assignment  of  its  structure  to  an  isolated  ion-pair. 


Experimental  Methods 

Helium  droplets  are  fonned  (10  ~  s'  )  by  the  continuous  expansion  of  He  gas  (30  bar) 
through  a  5  pm  diameter  pin-hole  nozzle.  The  nozzle  is  cooled  by  a  closed-cycle  refrigerator  to 
«14  K.  The  average  droplet  size  under  these  source  conditions  is  «6000  He  atoms.65, 66  Upon 

67  68 

leaving  the  high  pressure  region  of  the  expansion,  the  droplets  cool  by  evaporation  to  0.4  K.  ’ 
The  droplet  expansion  is  skimmed  into  a  beam,  which  then  passes  into  a  differentially  pumped 
“pick-up”  chamber.  The  droplets  are  doped  with  [emimjITfrN']  as  the  beam  passes  directly 
above  the  opening  to  a  resistively  heated  quartz  oven  containing  the  liquid.  Zaitsau,  et.  al.  used 
the  integral  effusion  Knudsen  method  to  obtain  [emim+][Tf2N‘]  vapor  pressures.30  Extrapolating 
their  data,  we  find  that  an  oven  temperature  of  440  K  produces  the  «1.5  x  10’5  Torr  vapor 
pressure  required  for  each  droplet  to  pick-up  one  [emim+][Tf2N']  ion-pair,  on  average.  The  IR 
spectra  are  measured  with  the  oven  source  at  420  K  to  reduce  the  vapor  pressure  slightly  below 
what  is  optimal  for  the  pick-up  of  a  single  ion-pair.  Under  these  conditions,  the  probability  of 
sequentially  picking  up  two  ion-pairs  into  the  same  droplet  is  less  than  10%  that  of  picking  up 
one,  and  the  consideration  of  the  cluster  spectra  can  be  neglected.  Previous  photoionization  mass 
spectrometry  24  and  photoelectron  spectroscopy  29  studies  of  [emim+][Tf2N']  used  an  effusive 
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oven  source  at  «470  K,  and  no  evidence  of  decomposition  was  found.  The  internal  energy  of 
He-solvated  [emim+][Tf2N‘]  is  removed  by  He  atom  evaporation  and  the  system  returns  to  0.4  K 
within  a  nanosecond  of  pick-up.69  During  this  process,  each  evaporating  He  atom  removes  on 
average  5  cm'1  (0.014  kcal/mol)  of  energy.69  Due  to  this  rapid  cooling  of  the  He  droplet  upon 
ion-pair  pick-up,  the  distribution  of  conformational  isomers  in  the  droplets  is  expected  to 
resemble  the  equilibrium  distribution  in  the  vapor,  as  has  been  observed  for  other  multi-isomeric 

*  57,63,64 

systems.  ’ 

The  He  droplets  are  detected  with  a  quadrupole  mass  spectrometer  equipped  with  an 
electron  impact  ionizer  (90  eV),  which  provides  signatures  of  both  the  neat  and  doped  fractions 
of  the  droplet  ensemble.  To  optimize  for  the  pick-up  of  single  ion-pairs,  the  oven  temperature  is 
adjusted  while  monitoring  the  mass  spectrum  of  the  droplet  beam,  as  described  below.  The  IR 
spectrum  of  the  He-solvated  [emim+][Tf2N'j  ion-pair  is  obtained  by  overlapping  the  droplet 
beam  with  the  output  from  a  continuous  wave  IR  Optical  Parametric  Oscillator  (OPO)  laser 
system.  This  overlap  is  achieved  by  either  counter-propagating  the  OPO  and  droplet  beams  or 
by  steering  the  OPO  beam  into  a  multipass  cell,  consisting  of  two  parallel  gold-coated  mirrors. 
Vibrational  excitation  of  the  solvated  [emim  jjTfxNr]  ion-pair  followed  by  vibrational  relaxation 
results  in  the  evaporation  of  -620  He  atoms  for  CH  stretch  excitation  near  3100  cm'1 .  This 
geometric  cross-section  reduction  leads  to  a  concomitant  decrease  in  the  electron  impact 
ionization  cross-section  of  the  [emim+][Tf2N‘]  doped  He  droplets  as  they  enter  into  the  mass 
spectrometer.  The  [emim  ][Tf2N"]  IR  spectrum  is  obtained  by  monitoring  the  depletion  of  the 
ion  signal  (nonnalized  to  laser  power)  as  the  OPO  counter-propagates  the  droplet  beam  and  is 
tuned  from  2800-3300  cm'1.  The  dipole  moment  measurements  described  below  are  obtained 
with  the  OPO  aligned  into  the  multipass  cell.  In  this  configuration,  the  entire  droplet  beam  / 
laser  interaction  region  is  contained  within  a  static  dc  electric  field  established  between  two 
parallel  stainless  steel  electrodes  situated  orthogonal  to  the  multipass  mirrors.  The  laser 
polarization  is  aligned  parallel  to  this  static  electric  field  with  a  Fresnel  Rhomb.  The  OPO 
wavelength  is  fixed  to  the  peak  of  an  IR  transition,  and  the  electric  field  is  linearly  ramped  from 
0-30  kV/cm.  The  field  ON  to  field  OFF  intensity  ratio  of  the  IR  band  is  averaged  for  about  10  of 
these  scans. 


6 

Distribution  A:  approved  for  public  release;  distribution  unlimited. 


Theoretical  Methods 


Previously,  Tsuzuki  and  co-workers  found  that  the  ab  initio  (MP2/6-31  1G**)  interaction 
energy  between  cation  and  anion  in  [emim+][Tf2N‘]  is  somewhat  insensitive  to  their  relative 
orientations,  and  they  found  16  conformational  isomers  with  electronic  energies  within  10 
kJ/mol  of  the  lowest  energy  species.42'43  We  have  carried  out  geometry  optimizations  and 
hannonic  frequency  computations  at  the  MP2/aug-cc-pVTZ  level71'75  of  theory  using 
GAMESS  ’  for  23  low  energy  isomers,  starting  from  the  structures  reported  in  Ref.  .  Here, 
the  harmonic  frequencies  are  scaled  by  0.96,  and  these  scaled  values  are  used  to  determine  zero- 
point  vibrational  energies  and  thennal  corrections.  The  five  lowest  energy  [emim  h][Tf2N']  ion- 
pair  isomers  at  this  level  are  shown  in  Fig.  1  and  numbered  1-5  according  to  their  relative  zero- 
point  corrected  electronic  energies.  The  relative  Gibbs  free  energies  computed  at  420  K  are  0.70, 
1.37,  0.00,  3.81  and  5.91  kJ/mol  for  isomers  1-5,  respectively.  The  entropic  effect  is  to  make 
isomer  3  the  most  abundant  configuration  at  thermal  equilibrium.  However,  we  emphasize  that 
the  relative  free  energies  of  these  five  isomers  are  within  the  error  bars  of  the  calculation  (±  «5 
kJ/mol),  and  it  is  therefore  difficult  to  know  with  any  certainty,  at  this  level  of  theory,  which 
isomer  is  actually  the  most  stable  one.  The  structures  for  the  complete  set  of  23  isomers  along 
with  their  relative  energetics  and  computed  spectra  are  given  in  the  Supporting  Information. 

The  vibrational  transition  moment  angles  (VTMAs)  ’  and  permanent  electric  dipole 
moment  vectors  (pp)  are  computed  at  this  same  level  of  theory  for  each  [emim  jfTfjN']  isomer. 
For  any  one  normal  mode  vibration,  the  VTMA  is  defined  as  the  angle  pp  makes  with  the 
transition  dipole  moment  vector  (pt).  Orientation  of  the  He-solvated  molecules  in  a  static  electric 
field  allows  for  the  simultaneous  measurement  of  pp  and  the  VTMA  associated  with  a  particular 
band  in  the  IR  spectrum.78, 79  This  is  achieved  by  measuring  the  field  dependence  of  the  IR  band 
intensity.  By  choosing  a  combination  of  VTMA  and  pp,  this  field  dependence  of  the  IR  intensity 
can  be  simulated  and  compared  to  the  experiment.  In  the  simulation,  the  average  projection  of  pt 
onto  the  OPO  polarization  axis  is  detennined  for  field  strengths  in  the  range  0-30  kV/cm.  The 
weighting  function  for  this  average  is  the  normalized  orientational  distribution  of  pp  about  the 
lab-frame  static  field  axis.  This  dipole  distribution  function,  P( cos0),  is  computed  for  each  static 
electric  field  strength  (increments  of  5  kV/cm)  according  to  the  procedures  of  Kong  and 
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Bulthuis.80' 81  Here,  0  is  defined  as  the  angle  between  pp  and  the  static  field  axis.  Calculations  of 
P{ cos0)  are  based  on  the  ab  initio  rotational  constants  (divided  by  3  to  account  for  the  helium 
solvent)57  and  inertial  dipole  moment  components  for  isomer  3  and  a  rotational  temperature  of 
0.4  K.  While  the  computed  orientational  distribution  is  sensitive  to  the  magnitude  of  pp,  a 
quantitative  detennination  of  rotational  constants  is  not  required  for  a  relatively  “heavy”  system, 
such  as  [emim  jfTfjN'].  Indeed,  the  simulation  will  be  insensitive  to  the  exact  values  of  the 
rotational  constants,  given  that  the  interaction  between  the  dipole  moment  and  the  electric  field  is 
much  larger  than  the  average  [emiml][Tf2N‘]  rotational  energy  at  0.4  K.  With  the  computed 
P( cos0)  distribution,  pt  is  projected  onto  the  polarization  axis  of  the  OPO  (parallel  to  the  lab- 
frame  field  axis)  and  averaged  over  0  to  give  the  relative  intensity,  Apara,  due  to  molecular 
orientation  at  a  given  field  strength,  ~  where  a  is  the  VTMA: 


Apara(a )  x  J  -f*(cos0  )[sin^ a  +  2cos20  -3sin2a  cos70]<i cos0  . 


(1) 


Results  and  Discussion 
A.  Mass  Spectrometry 

Electron  impact  ionization  of  a  He  droplet  produces  a  solvated  He+  cation,  which  leads  to 
a  vibrationally  excited  He2+  cation.  The  outcome  of  this  ionization  process  is  now  well  known  to 
produce  a  distribution  of  gas  phase  Hen+  ions.83  The  mass  spectrum  in  Fig.  2a  shows  the  Hen+ 
ions  associated  with  the  ionization  of  the  neat  droplet  beam.  If  the  droplets  are  doped,  in 
addition  to  the  Hen+  distribution,  ions  associated  with  the  solvated  molecule  can  be  produced  as  a 
result  of  the  charge  transfer  reaction, 

He+  +  dopant  — >  He  +  [dopant]+  — »  He  +  [dopant  fragments]4)  (2) 

Figures  2b-d  show  the  evolution  of  the  mass  spectrum  as  the  quartz  oven  source  is  heated  to 
dope  the  droplets  with  [emim  jfTfjlSr].  Due  to  the  ionization  of  what  is  presumably  the  He- 
solvated  [emim  ][Tf2N‘]  ion-pair,  heating  the  oven  to  408  K  results  in  a  new  peak  in  the  mass 
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spectrum  at  m/z=l  1 1 .  This  peak  corresponds  to  the  non- fragmented  [emim  ]  cation.  Upon 
heating  the  oven  further  to  426  K,  the  peak  at  m/z=T  1 1  becomes  more  intense  and  other  peaks  in 
the  lower  mass  range  begin  to  appear.  The  neat  He  droplet  mass  spectrum  (Fig.  2a)  is  subtracted 
from  the  mass  spectrum  recorded  with  the  oven  at  426  K  (Fig.  2d),  giving  the  difference 
spectrum  shown  in  Fig.  3.  With  the  He  peaks  subtracted,  the  remaining  peaks  are  due  only  to 
ions  formed  upon  He+  charge  transfer  to  species  picked-up  by  the  droplet  beam.  In  addition  to 
m/z=l  11,  other  predominant  peaks  in  the  mass  spectrum  such  as  m/z=3 1,  50,  64,  69  are  due  to 
CF+,  CF2+,  SCW,  and  CF3+,  respectively,  which  are  fragments  from  [Tf2N'].  In  previous  line  of 
sight  electron  impact  ionization  mass  spectrometry  and  VUV  soft  photoionization  studies," 
the  dominance  of  the  m/z=l  1 1  [emim  ]  peak  was  interpreted  as  being  due  to  the  dissociative 
ionization  of  [emim  ][Tf2N‘]  ion-pairs  in  the  gas  phase,  according  to  the  mechanism 

[emim+]  [Tf2N‘]  +  hv  — »  [emim+]  [Tf2N]  +  e'  — >  [emim+]  +  [Tf2N]  +  e'.  (3) 

A  similar  dissociative  ionization  mechanism  is  likely  operative  upon  electron  impact  ionization 
of  [emim][Tf2N']  doped  He  droplets.  This  appears  reasonable  because  the  ionization  of 
[emim  ][Tf2N']  is  only  likely  to  proceed  via  He+  charge  transfer,  and  the  difference  in  the  He  and 
[emim  ][Tf2N']  ionization  potentials  is  «15  eV,  which  is  much  larger  than  the  «0.5  eV  binding 
energy  of  the  cation-radical  complex  formed  upon  ionization.  Therefore,  the  dominance  of 
m/z=l  1 1  in  the  He  droplet  mass  spectrum  provides  strong  evidence  for  the  efficient  pick-up  of 
[emim+][Tf2N‘]  ion-pairs. 


B.  Infrared  Spectroscopy 

The  IR  spectrum  shown  in  Figure  4  (middle  trace)  is  obtained  by  monitoring  the  ion 
depletion  signal  in  mass  channel  m/z=l  1 1  as  the  OPO  is  continuously  tuned  from  2800  to  3250 
cm'1.  There  are  several  spectral  features  in  this  wavelength  region,  presumably  due  to  the  CH 
stretches  of  [emim1  ][Tf2N']  ion-pairs.  The  ~2  cm'1  linewidth  of  each  band  is  broader  than  the 
expected  rotational  contour  at  0.4  K,  and  this  width  is  therefore  attributed  to  the  vibrational 
relaxation  timescale  of  the  excited  molecule  (i.e.  «3  ps).  The  red  trace  at  the  top  of  Fig.  4  is  the 
Boltzmann  weighted  sum  of  the  ab  initio  hannonic  spectra  (weights  determined  from  free 
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energies  at  420  K).  There  are  three  distinctive  CH  stretching  regions  of  the  computed  harmonic 
spectra,  and  these  correspond  to  the  symmetric  stretching  motions  of  the  imidazolium  (IM)  alkyl 
side  chains  (2950  to  3000  cm'1),  the  antisymmetric  stretching  motions  of  the  alkyl  side  chains 
(3000  to  3 100  cm'1),  and  the  three  CH  stretches  on  the  IM  ring  (3090  to  3200  cm'1).  The  IM  ring 
CH  stretch  modes  consist  of  the  symmetric  and  antisymmetric  stretching  motions  of  the  H-C(4)- 
C(5)-H  group  («3190  cm'1  and  ~3175  cm"1,  respectively)  and  the  C(2)-H  stretch.  The  computed 
CH  stretches  are  rather  weak,  with  the  majority  having  intensities  less  than  15  km/mol.  The 
exceptions  to  this  are  the  C<2)-H  stretching  bands  associated  with  those  isomers  that  have 
C(2)H+  anion  H-bonded  configurations  (e.g.  isomer  3).  For  these  configurations,  the  strong,  in¬ 
plane  H-bonds  lead  to  characteristic  red  shifts  and  intensity  enhancements  for  the  C(2)-H 
stretching  bands.  The  predicted  C(2)-H  band  origins  range  from  «3090  cm'1  for  the  most  strongly 
H-bonded  isomer  (isomer  3)  to  «3160  cm"1  for  the  “non-H-bonded”  isomers  (isomers  1  and  2). 
The  C(2)-H  stretching  bands  are  labeled  with  bold  numbers  above  the  red  trace  in  Fig.  4  to 
indicate  the  associated  isomer. 

Although  the  weighted  ab  initio  spectrum  allows  for  a  qualitative  assignment  of  the 
experiment,  there  is  a  degree  of  spectral  complexity  in  the  H-bonded  C(2)-H  stretch  region  that  is 
not  reproduced  in  the  harmonic  frequency  computations.  Indeed,  the  number  of  transitions  and 
their  relative  positions  and  intensities  are  not  reproduced  in  the  3050-3150  cm'1  region.  One 
possible  origin  of  this  discrepancy  is  the  accuracy  associated  with  the  relative  free  energies  used 
to  weight  the  simulated  spectrum.  It  is  not  obvious  a  priori  whether  or  not  the  calculations  are 
capable  of  determining  the  relative  energetics  for  these  lowest  five  isomers  to  an  accuracy  of  a 
few  kJ/mol.  Another  potential  source  of  error  is  the  harmonic  approximation  employed  to 
compute  the  vibrational  spectrum.  Lassegues  and  co-workers  computed  (B3LYP/6-31G**) 
cubic  force  constants  for  the  bare  [emim+]  cation  and  suggested  that  the  relatively  large 
magnitudes  for  some  of  these  should  result  in  anharmonic  resonances  in  the  CH  stretch  region 
for  IM  based  RTILs.  Indeed,  these  computations  do  indicate  that  the  IM  CH  stretches  are 
coupled  to  the  overtones  and  combination  of  two  IM  ring  vibrations  near  1600  cm"1,  as  the  cubic 
force  constants  associated  with  these  couplings  are  particularly  large  (~70  cm"1)  and  the  energy 
differences  of  the  modes  involved  are  relatively  small.  This  motivates  an  analysis  of  the 
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observed  vibrational  complexity  in  the  3050-3150  cm"1  region  in  terms  of  anhannonic 
resonances  between  the  IM  ring  stretching  overtones/combination  and  the  IM  CH  stretches. 


The  spectral  features  associated  with  these  predicted  anhannonic  resonances  can  be 
estimated  with  a  reduced  dimensional  model  that  treats  the  anhannonic  couplings  explicitly  by 
diagonalizing  the  6x6  coupling  matrix  shown  in  eq.  4.  Here,  Vi  through  V3  are  the  fundamental 
symmetric  C(4>5)-H,  antisymmetric  C^^-H,  and  C(2)-H  stretching  modes,  respectively;  v4  and  v5 
are  the  fundamental  antisymmetric  and  symmetric  ring  breathing  modes,  respectively.  The 
diagonal  elements  are  the  harmonic  frequencies  (MP2/aug-cc-pVTZ),  which  have  been  scaled  by 
0.96  to  account  for  diagonal  anhannonic ity.  In  this  model,  the  potential  expansion  is  truncated  at 
the  cubic  terms,  and  the  only  non-negligible  interactions  are  those  that  couple  the  high  frequency 
CH  stretch  fundamentals  to  the  overtones/combination  of  the  two  lower  frequency  ring  modes. 
Moreover,  any  interactions  involving  the  alkyl  CH  stretch  modes  are  neglected,  because  these 
are  expected  to  be  only  weakly  coupled  to  the  vibrations  in  the  ring.  The  cubic  force  constants 
are  taken  from  Ref.  ,  and  the  multiplicative  factors  in  the  off-diagonal  matrix  elements 
originate  from  the  matrix  elements  of  the  nonnal  coordinates  in  the  harmonic  oscillator  basis.  ’ 

85 

This  model  further  assumes  that  the  cubic  force  constants  do  not  depend  on  the  presence  of 
the  anion,  given  that  they  were  computed  only  for  the  bare  [emim  ]  cation.'  The  intensities  of 
the  fundamental  CH  stretch  basis  modes  are  those  computed  with  the  harmonic  approximation, 
and  the  zeroth-order  intensities  for  2v4,  2vs  and  V4+V5  are  set  to  zero,  which  is  justified  given  the 
hannonic  nature  of  the  potential  energy  surface  along  the  v4  and  V5  nonnal  coordinates.50  For 
simplicity,  we  estimate  the  anhannonic  intensities  as  ,  in  which  i  is  the  harmonic 

intensity  vector  and  U  is  a  matrix  containing  the  squares  of  the  eigenvector  coefficients.  This 
treatment  is  carried  out  for  isomers  1-5,  and  the  resulting  spectra  are  Boltzmann  weighted  and 
summed  to  give  the  blue  simulation  shown  at  the  bottom  of  Fig  4.  In  this  anhannonic 
simulation,  the  transition  frequencies  and  intensities  for  the  alkyl  group  CH  stretches  are  the 
same  as  those  in  the  harmonic  simulation.  The  individual  anhannonic  spectra  for  each  isomer 
are  given  in  the  Supporting  Infonnation. 

This  simple  approximate  model  leads  to  a  noticeable  improvement  in  the  agreement 
between  the  experiment  and  the  number  of  intense  bands  predicted  in  the  simulation,  which 
again  assumes  a  thermal  distribution  of  ion-pair  isomers.  The  CH  stretch  transitions  above  3150 
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cm’1  are  similar  to  those  in  the  harmonic  simulation,  with  the  anharmonic  coupling  resulting  in  a 
5-10  cm’1  blue  shift  for  each  of  these  bands.  This  relatively  modest  perturbation  of  the  C^^-H 
stretch  vibrations  is  due  to  the  large  energy  differences  between  these  modes  and  the  ring 
overtones/combination.  Similarly,  the  non-H-bonded  C(2)-H  stretches  of  the  “stacked”  isomers 
(1,  2  and  5)  are  largely  unperturbed.  One  or  more  these  isomers  are  likely  present  in  the 
distribution,  because  the  broad  feature  near  3150  cm’1  agrees  well  with  the  predicted  C(2)-H 
stretch  origins  for  these  species.  In  contrast  to  these  higher  frequency  modes,  significant 
perturbations  are  observed  for  the  lower  frequency  C(2)-H  stretch  modes  associated  with  isomers 
containing  strong,  directional  C(2)H+  anion  H-bonds.  The  H-bonding  interaction  shifts  the  C(2)- 
H  stretch  to  lower  energy  and  into  resonance  with  the  ring  overtones/combination,  resulting  in 
significant  mode  mixing  and  intensity  borrowing.  For  example,  the  predicted  C(2)-H  stretch 
bands  for  isomers  3  and  4  (scaled  harmonic  3093  and  3107  cm’1,  respectively)  are  particularly 
intense  due  to  H-bonding,  and  the  interaction  between  these  modes  and  the  ring 
overtone/combination  states  results  in  a  characteristic  four  peak  pattern  for  each  isomer.  This 
predicted  pattern  for  these  H-bonded  isomers  is  consistent  with  that  observed  experimentally 
between  3050  and  3125  cm'1.  We  emphasize,  however,  that  the  band  positions  and  intensities 
are  rather  sensitive  to  both  the  magnitude  of  the  coupling  constants  and  the  zeroth-order  band 
origins,  and  we  therefore  do  not  expect  to  achieve  quantitative  agreement  between  the 
experiment  and  simulation  at  this  level  of  theory  and  with  this  approximate  model.  The  arrows 
and  brackets  in  Fig.  4  indicate  the  suggested  assignments.  In  general,  the  agreement  between  the 
Boltzmann  weighted  anharmonic  simulation  and  the  experimental  IR  spectrum  suggests  that  the 
droplets  pick-up  and  cool  a  gas-phase  distribution  of  ion-pair  isomers  (both  H-bonded  and  non- 
H-bonded),  which  are  at  thennal  equilibrium  at  the  «420  K  source  temperature;  otherwise,  at  0.4 
K,  only  the  lowest  free  energy  isomer  would  be  expected. 
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C.  Polarization  Spectroscopy 

The  electric  field  dependence  of  the  band  centered  at  3 1 13  cm’1  (indicated  by  *  in  Fig.  4) 
is  measured  as  the  ratio  of  the  field  ON  to  field  OFF  signal  over  a  range  of  applied  field 
strengths,  and  this  is  shown  as  the  black  trace  in  Fig.  5.  This  band  was  selected  because  it  does 
not  significantly  overlap  with  other  bands,  and  we  assign  this  feature  to  the  C(2)-H  stretch 
transition  of  isomer  3.  For  isomer  3,  the  squares  of  the  eigenvector  coefficients  for  the  upper 
level  of  this  transition  are  c2vi=0.006,  c2V2=0.003,  c2V3=0.54,  c22V4=0.18,  c22V5=0.037,  and 
c'v4+v5=0.23,  indicating  the  substantial  mode  mixing  described  above.  The  field  ON:OFF 
intensity  ratio  is  simulated  by  assuming  the  ab  initio  pp  for  isomer  3  (10.8  D),  and  this  simulation 
is  shown  as  the  red  squares  in  Fig.  5.  Also  shown  are  the  simulations  for  three  hypothetical 
structures  having  dipole  moments  equal  to  5,  9  and  15  D.  For  all  four  simulations,  the  VTMA 
(a)  is  set  to  37°,  which  is  chosen  to  reproduce  the  asymptotic  (high-field  limit)  field  dependence. 
We  note  that  because  the  3113  cm'1  band  is  a  member  of  an  anhannonic  resonance  polyad,  the 
experimentally  derived  VTMA  (37°)  cannot  in  general  be  compared  to  the  ab  initio  VTMA 
obtained  from  the  harmonic  frequency  computation  («27  0  for  the  C(2)-H  stretch  of  isomer  3),  and 
the  VTMA  here  serves  only  as  a  single  adjustable  fitting  parameter.  In  the  limit  of  infinite  field, 
the  asymptotic  field  dependence  is  due  only  to  the  VTMA  (a),  and  is  independent  of  pp  (as  long 


13 

Distribution  A:  approved  for  public  release;  distribution  unlimited. 


as  pp  is  finite).  For  example,  when  a=0°,  54.75°  or  90°  (parallel  polarization)  the  high  field  limit 
of  the  ON:OFF  intensity  ratios  are  3,  1,  or  0,  respectively.  However,  the  rate  at  which  these 
asymptotic  limits  are  approached  depends  strongly  on  the  magnitude  of  pp.  Indeed,  the  pp=5  D 
system  approaches  more  slowly  the  asymptotic  limit  than  that  with  a  pp  equal  to  9  D.  Clearly, 
there  is  excellent  agreement  between  the  experimental  trace  and  the  simulation  that  assumes  an 
isomer  3  structure  and  permanent  dipole  moment.  From  this  polarization  spectroscopy 
measurement,  we  assign  a  dipole  moment  of  1 1±2  D  to  the  carrier  of  this  IR  band.  Again,  this 
agrees  quite  well  with  the  predicted  pp  for  isomer  3,  and  more  generally,  it  agrees  with  the  range 
of  dipole  moments  (9-13  D)  that  are  predicted  for  all  [emim4][Tf2N"]  ion-pair  isomers.  The 
computed  magnitudes  of  pp  for  isomers  1-5  are  9.9,  9.7,  10.8,  10.7  and  1 1.0,  respectively. 

Indeed,  this  result  provides  rather  definitive  confirmation  that  ion-pairs  are  the  main  constituents 
of  [emim+][Tf2N']  vapor  at  420  K. 

Previous  VUV  photoionization  studies  of  [emim+][Br‘]  showed  that  vaporization  occurs 
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via  the  evolution  of  alkyl  bromides  and  alkylimidazoles,  rather  than  the  intact  ion-pairs. 
Experimental  activation  enthalpies  were  found  to  be  in  good  agreement  with  computations  of  the 

formation  of  these  species,  in  which  the  ion-pair  thennal  decomposition  pathway  consists  of  an 

86 

Sn2  reaction  involving  the  bromide  anion  and  either  of  the  alkyl  substituents  on  the  IM  ring. 

The  decomposition  of  [emim+][Tf2N']  via  this  mechanism  would  produce  methylimidazole  or 
ethylimidazole,  which  have  computed  dipole  moments  of  4.2  and  4.3  D,  respectively.  Another 
feasible  decomposition  pathway  is  proton  transfer  from  [emim]  to  [Tf2N‘]  to  produce  l-ethyl-3- 
methylimidazole  2-carbene  (2.4  D  dipole  moment).  ’  In  general,  decomposition  products  are 
expected  to  have  dipole  moments  that  are  qualitatively  smaller  than  those  computed  for  the  intact 
ion-pair.  The  polarization  spectroscopy  reported  here  indicates  that  these  decomposition 
pathways  for  [emim  ][Tf2N']  are  negligible  at  420  K,  and  its  vaporization  at  this  temperature 
largely  produces  intact  ion-pairs. 
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Conclusions 


The  ionic  liquid  1 -ethyl-3 -methylimidazolium  bis(trifluoromethylsulfonyl)imide  was 
vaporized  at  420  K,  solvated  in  0.4  K  He  nanodroplets  and  probed  with  IR  laser  spectroscopy. 
The  experimental  conditions  were  such  that  each  He  droplet  picked  up  a  single  molecule  from 
the  gas  phase.  The  IR  spectrum  in  the  CH  stretch  region  was  measured,  and  it  reveals  bands  that 
are  assigned  to  intact  [emim+][Tf2N‘]  ion-pairs.  The  assignment  is  based  on  a  comparison  of  the 
IR  spectrum  to  ab  initio  harmonic  frequency  computations  of  23  ion-pair  isomers  that  were 
found  to  be  local  minima  at  the  MP2/aug-cc-pVTZ  level  of  theory. 

The  agreement  between  the  experimental  and  computed  IR  spectra  is  improved  by  taking 
into  account  the  anhannonic  coupling  between  the  IM  CH  stretch  fundamentals  and  the 
overtones  and  combination  of  two  IM  ring  stretching  vibrations  near  1600  cm'1.  For  isomers 
having  a  strong  H-bond  between  the  cation  and  anion,  the  H-bonded  donor  C(2)H  stretch  on  the 
IM  ring  is  shifted  to  lower  energy  and  into  resonance  with  the  ring  stretching  overtones  and 
combination.  The  anhannonic  coupling  between  these  modes  results  in  a  resonance  polyad 
consisting  of  four  peaks  centered  near  3075  cm'1  that  are  spread  over  «100  cm"1.  The 
experimental  IR  spectrum  agrees  qualitatively  with  this  characteristic  four  peak  pattern  predicted 
for  the  H-bonded  ion-pair  isomers.  This  pattern  of  transitions  is  absent  for  isomers  lacking  a 

strong,  directional  C(2)H+ "anion  H-bond,  and  the  C(2)H  stretching  band  in  these  “stacked” 
isomers  are  predicted  as  single  transitions  in  the  3 125-3175  cm'1  range.  Vibrational  bands  are 
observed  in  this  wavelength  range,  indicating  that  non-H-bonded,  stacked  isomers  are  also 
present  in  the  ensemble  of  droplets.  The  experimental  IR  spectrum  is  therefore  consistent  with 
the  pick-up  and  cooling  of  ion-pairs  by  He  droplets,  in  which  the  cooling  is  sufficiently  fast  to 
freeze  out  the  thermal  distribution  of  the  vapor  (at  420  K). 

The  assignment  to  intact  [emim  jfTf^N']  ion-pairs  is  confirmed  with  polarization 
spectroscopy,  in  which  the  permanent  electric  dipole  moment  of  the  He-solvated  species  is 
oriented  in  a  lab-frame  dc  electric  field.  The  electric  field  dependence  of  several  IR  band 
intensities  is  consistent  with  a  He-solvated  molecule  that  has  a  dipole  moment  in  the  range  1 1+2 
D.  The  computed  magnitudes  of  qp  for  [emim+][Tf2N']  ion-pair  isomers  are  indeed  between  9 
and  13  D,  and  the  lowest  free  energy  isomer  computed  at  the  MP2/aug-cc-pVTZ  level  has  a  10.8 
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D  dipole  moment.  The  potential  products  of  thennal  decomposition  have  computed  dipole 
moments  that  are  considerably  lower  (<4.3  D)  than  those  found  for  the  intact  ion-pairs.  The 
polarization  spectroscopy  employed  here  therefore  indicates  that  the  thennal  decomposition  of 
[emim  ][Tf2N‘]  is  negligible  at  420  K,  and  that  its  vaporization  at  this  temperature  largely 
produces  intact  ion-pairs. 

The  results  presented  here  are  inconsistent  with  previous  interpretations  of  the  IR  spectra 
of  liquid  [emim+][Tf2N'],  in  which  the  temperature,  pressure  and  dilution  effects  on  the  hydrogen 
bonding  network  in  the  liquid  were  analyzed  in  terms  of  the  number  and  relative  intensities  of 
bands  in  the  IM  CH  stretch  region,  without  consideration  of  anhannonic  polyads.46  As  originally 
suggested  by  Lassegues  and  co-workers,'  anhannonic  resonances  in  the  CH  stretch  region 
cannot  be  neglected  in  these  analyses.  However,  contrary  to  previous  suggestions,  H-bonding 
interactions  in  the  liquid  are  also  essential  to  explain  the  complexity  in  the  3050-3150  cm’1 
region  of  the  condensed  phase  IR  spectra,  as  discussed  by  Ludwig  and  co-workers.50  Indeed,  this 
is  now  particularly  evident  from  the  spectra  of  the  cold,  He  nanodroplet  isolated  species,  in 
which  both  H-bonding  interactions  and  anhannonic  coupling  to  overtones  and  combination  states 
must  be  invoked  to  rationalize  the  spectra. 
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Figure  Captions 


Figure  1:  Five  low  energy  [emim+][Tf2N’]  ion-pair  structures  (electronic  +ZPVE)  at  the 
MP2/aug-cc-pVTZ  level  of  theory. 

Figure  2:  a)  Mass  spectrum  of  the  He  droplet  beam  with  the  [emim+][Tf2N’]  source  off.  Peaks  at 
multiples  of  4  u  are  due  to  Hen+  cluster  ions.  The  peak  at  m/z=18  u  is  due  to  a  small  fraction  of 
the  droplets  (<1%)  that  have  picked  up  residual  water  molecules  in  the  vacuum  chamber.  The 
mass  spectrum  of  the  droplet  beam  with  the  [emim+][Tf2N']  source  heated  to  b)  408  K,  c)  418  K, 
and  d)  426  K  shows  the  gradual  appearance  of  the  intact  [emim  ]  cation  at  m/z=l  1 1  u. 
Temperatures  are  obtained  with  a  K-type  thermocouple  attached  to  the  quartz  furnace. 

Figure  3:  The  mass  spectrum  of  the  neat  droplet  beam  was  scaled  and  subtracted  from  the  mass 
spectrum  with  426  K  oven  temperature  conditions  in  order  to  reduce  contributions  from  the  Hen+ 
peaks.  In  addition  to  the  [emim+]  peak  at  m/z=l  1 1  u,  peaks  at  3 1,  50,  64,  and  69  u  are  observed 
and  are  due  to  CF+,  CF2+,  S02+,  and  CF3+,  which  are  fragments  from  [Tf2N'].  Peaks  at  m/z=27, 
28,  and  29  u  are  characteristic  fragment  peaks  of  [emim1]  associated  with  the  ethyl  side  chain. 

Figure  4:  Helium  Droplet  IR  spectrum  in  the  CH  stretch  region  for  [emim  ][Tf2N’]  vapor.  The 
He  droplet  spectrum  corresponds  to  the  wavelength  dependent  ion  signal  depletion  in  the 
m/z=T  1 1  u  channel.  The  top  (red)  simulated  spectrum  is  a  Boltzmann  weighted  sum  of  the 
computed  harmonic  frequency  spectra  (scaled  by  0.96)  for  the  five  lowest  free  energy  isomers 
computed  at  the  MP2/aug-cc-pVTZ  level.  The  bottom  (blue)  simulated  spectrum  is  a  Boltzmann 
weighted  sum  of  the  computed  anhannonic  frequency  spectra  for  these  same  five  isomers.  The 
peaks  marked  by  #  are  the  four  transitions  associated  with  the  anhannonic  resonance  polyad  that 
involves  the  H-bonded  C(2)-H  stretch  and  the  IM  ring  stretching  overtones  and  combination. 

Figure  5:  Field  ON:OFF  intensity  ratio  for  the  3113  cm’1  band  as  a  function  of  the  applied 
orientation  field.  The  laser  polarization  is  aligned  parallel  to  the  static  field,  and  the  asymptotic 
dependence  of  the  signal  indicates  a  VTMA  close  to  «37  degrees.  The  smooth  lines  are 
simulations  of  the  field  dependence  for  several  different  magnitudes  of  the  permanent  electric 
dipole  moment,  as  given  in  the  inset.  The  rate  at  which  the  signal  approaches  the  asymptotic 
limit  is  consistent  with  a  pennanent  electric  dipole  moment  of  1 1±2  D. 
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